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ABSTRACT 

This  paper  shows  how  the  principle  of  neutron  radiog¬ 
raphy  can  be  used  to  indirectly  quantify  the  critical  physi¬ 
cal  state  of  lithium  concentration  across  battery  electrodes 
at  steady-state  conditions  (after  a  long  relaxation  time  or 
small  load)  as  a  first  step  in  this  important  effort  to  measure 
in-situ  battery  physical  states  and  validate  electrochemical 
battery  models.  A  model  of  the  expected  loss  in  beam  in¬ 
tensity  after  passing  through  the  different  layers  of  a  battery 
pouch  cell  is  constructed  based  on  the  material  densities 
and  dimensions.  This  model  is  augmented  with  simulation 
of  the  neutron  transmission  behavior,  including  optical  ef¬ 
fects  due  to  the  geometric  unsharpness  and  the  detector  re¬ 
sponse.  The  resulting  model  provides  the  basis  for  a  com¬ 
prehensive  simulation  of  the  in-situ  metrology  of  lithium 
concentration  in  Li-ion  batteries,  and  comparison  with  ex¬ 
perimental  results. 

1.  INTRODUCTION 

Battery  management  requires  accurate  prediction  of 
bulk  and  the  spatiotemporal  Lithium  (Li)  concentrations. 
The  bulk  concentration  is  a  good  metric  for  the  total  avail¬ 
able  energy  stored  in  the  battery.  The  spatial  distribu¬ 
tion  is  important  for  predicting  the  available  power.  The 
minimum  and  maximum  local  concentrations  of  Li-ions 
throughout  the  electrode  and  the  electrolyte  solution  limit 
the  performance  and  longevity  of  the  battery  (Smith,  2010; 
Domenico  et  ak,  2010).  The  electrode  material  can  be 
modeled  as  a  distribution  of  spherical  particles  along  the 
electrode,  as  shown  in  Fig.  1.  The  ratio  of  Li  surface 
concentration  over  the  maximum  concentration  represents 
the  amount  of  readily  available  power  because  the  surface 
lithium  is  easily  removed  (Santhanagopalan  et  ak,  2008; 
Speltino  et  ak,  2009;  Smith  et  ak,  2010),  whereas  the  ions 
near  the  center  of  the  particle  are  limited  by  the  rate  of  dif¬ 
fusion  to  reach  the  surface. 

Hence,  our  ability  to  control  the  power  and  fully  utilize 
the  energy  from  Li-ion  batteries  depends  on  physics-based 


Figure  1:  Pouch  cell  construction  showing  an  expanded 
view  of  the  electrode  structure.  The  Lithium  concentration 
distribution  varies  across  the  electrode  and  along  the  radius 
of  the  spherical  particles  which  are  disturbed  across  the 
electrode,  during  battery  charging.  Modified  from  (Smith, 
2010;  Speltino  et  ak,  2009). 


and  scalable  battery  models  of  the  spatiotemporal  varia¬ 
tions  in  lithium  concentration  throughout  the  electrode  and 
electrolyte.  Although  there  are  many  distributed  parame¬ 
ter  micro-macro  scale  physics  based  models  of  the  elec¬ 
trochemical  reactions  and  mass  transport  phenomena  oc¬ 
curring  inside  the  battery,  they  could  all  benefit  from  mea- 
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surement  of  the  lithium  concentrations  for  model  valida¬ 
tion.  Neutron  radiography  offers  a  unique  opportunity  for 
measuring  in-situ,  the  changes  in  lithium  concentration  and 
hence  can  be  used  for  parameterization  and  validation  of 
mathematical  models  of  Li-ion  cells.  Prior  efforts  to  utilize 
neutron  imaging  resulted  in  only  qualitative  results  (Ka- 
mata  et  al.,  1997;  Lanz  et  al.,  2001). 

This  manuscript  includes  data  and  a  metrology  simu¬ 
lation  that  fill  the  critical  gap  of  experimental  verification 
and  availability  of  in-situ  lithium  concentration  patterns  in 
the  area  of  lithium-ion  batteries.  The  technological  break¬ 
through  from  the  availability  of  validated  models  will  make 
a  substantial  impact  to  the  hybrid  system  level  sizing  and 
power  management  for  an  ultra-light  far-reaching  portable 
power  source.  In  this  paper,  we  investigate  the  capabilities 
of  neutron  imaging  for  in  situ  measurement  of  the  spatial- 
temporal  concentration  of  Li  in  a  battery.  We  discuss  the 
physics  of  the  image  formation  process  using  pinhole  op¬ 
tics,  the  Lambert-Beer  law,  and  the  energy  dependent  neu¬ 
tron  cross  sections.  The  pouch  cell  construction,  shown  in 
Fig.  1,  consisting  of  multiple  layers  was  imaged  using  the 
thermal  neutron  imaging  beamline  at  the  National  Institute 
for  Standards  and  Technology  (NIST).  Experiments  were 
conducted  with  the  high  resolution  Micro-Channel  Plate 
(MCP)  detector,  which  is  a  state  of  the  art  neutron  detec¬ 
tor  at  the  NIST  Center  for  Neutron  Research  (NCNR). 

The  experiments  were  carefully  designed  to  achieve 
high  measurement  confidence.  Measurement  uncertainty 
arises  primarily  from  the  random  nature  of  the  neutron 
emission  and  interactions,  where  the  number  of  neutrons 
counted  in  a  given  period  of  time  can  be  described  by  a 
Poisson  distribution.  In  order  to  minimize  the  measurement 
uncertainty,  we  should  increase  the  number  of  counted  neu¬ 
trons,  which  can  be  achieved  by  careful  choice  of  materials, 
increasing  the  exposure  time,  and/or  by  spatial  averaging. 


utilization  of  neutron  imaging  for  Fuel  Cell  (FC)  liquid  wa¬ 
ter  metrology  has  been  demonstrated  at  imaging  facilities 
by  NIST  (Hussey  et  al.,  2010)  and  PSI  (Lehmann  et  al., 
2009),  and  advanced  our  understanding  of  the  water  dy¬ 
namics  and  associated  FC  degradation  phenomena  (Hick- 
ner  et  al.,  2008;  Siegel  et  al.,  2008).  A  similar  approach  can 
be  applied  to  measure  the  Li  movement  between  the  pos¬ 
itive  and  negative  electrodes  of  an  operating  Lithium  ion 
battery,  since  Li  also  has  a  large  cross  section  for  thermal 
neutrons.  Specifically  the  areas  shown  in  Fig.  2  indicate 
the  probability  of  neutron  interaction  with  atoms  of  various 
materials  used  in  the  construction  of  Lithium-ion  Batteries. 
When  a  neutron  interacts,  it  has  a  much  smaller  probability 
of  reaching  the  detector,  which  leads  to  beam  attenuation. 
The  change  in  Li  concentration  causes  a  change  in  neu¬ 
tron  beam  intensity  which  reaches  the  detector  after  passing 
through  the  battery,  as  shown  in  Fig.  3.  The  largest  change 
in  molar  concentration  of  lithium  occurs  within  the  inter¬ 
calation  materials.  Hence,  the  change  in  detected  beam  in¬ 
tensity  can  be  related  to  the  change  in  Li  concentration  in 
the  solid  electrodes. 


Figure  2:  Averaged  total  (scattering  and  absorbing)  neutron 
cross-sections  based  on  naturally  occurring  abundance  of 
isotopes. 


The  goal  of  this  paper  is  to  demonstrate  the  obser¬ 
vation  and  quantification  of  lithium  concentration  at  var¬ 
ious  states  of  charge.  Our  long-term  goal  is  to  observe 
the  spatiotemporal  distribution  of  lithium  across  the  two 
electrodes  during  charging  and  discharging  at  various  C- 
rates  (or  current  densities)  using  the  application  of  peri¬ 
odic  charge-discharge  profiles  and  stroboscopic  averaging 
to  overcome  the  spatio-temporal  resolution  limitation  of  the 
imaging  equipment.  In  this  paper  we  show  the  feasibility  of 
in-situ  Lithium  concentration  measurements  along  the  elec¬ 
trode  thickness  and  length,  at  various  states  of  charge  and 
discharge  by  imaging  steady-state  conditions,  i.e.  when  the 
battery  is  at  rest. 

2.  LITHIUM  METROLOGY 

Neutron  imaging  is  a  valuable  tool  for  in-situ  measure¬ 
ments,  similar  to  x-ray  radiography.  The  change  in  atten¬ 
uation  of  a  neutron  beam  passing  through  a  sample  can  be 
used  to  measure  a  change  in  the  material  over  time.  The 


2.1  Principles  of  Neutron  Imaging 

The  attenuation  of  neutron  beam  intensity  is  described 
by  the  Lambert-Beer  Law 

/  =  /0exp(-jic5)  (1) 

where  Iq  is  the  incident  intensity  which  passes  through  a 
material  of  width  8  (cm)  with  attenuation  coefficient  /i. 
Neutron  interactions  can  occur  as  either  scattering  or  ab¬ 
sorbing  or  both,  depending  on  the  material.  For  the  pur¬ 
pose  of  imaging  we  consider  the  effects  of  both  scatter¬ 
ing  and  absorption  since  both  effects  will  prevent  neutrons 
from  reaching  the  detector,  by  simply  adding  the  scatter¬ 
ing  and  absorbing  cross  sections  for  each  material  to  get  an 
estimate  of  the  total  beam  attenuation. 

As  a  first  order  approximation  for  an  object  consisting 
of  many  materials  the  total  attenuation  can  be  calculated 
from  the  composite  neutron  interaction  probably  of  each 
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independent  material  scaled  with  the  molar  concentration 
of  that  material.  Since  the  battery  is  composed  of  many 
materials,  we  can  write  the  total  linear  attenuation  as  the 
sum  of  two  components,  one  which  is  changing  over  time, 
the  lithium  concentration,  cu{t)  mol  cm-3,  and  one  which 
remains  stationary; 

I{t)  =Iq exp ( — OijCm (t )NA 8  —  L  a,c,A fAS)  (2) 

i,  i^Li 

where  c  is  the  molar  concentration  (mol  cm-3),  NA  is  Avo- 
gadro’s  Number,  and  ay,  is  the  neutron  cross  section  (cm2), 
shown  in  Fig.  2.  The  tabulated  neutron  cross  sections  for 
each  element  (isotope)  can  be  found  on  the  NCNR  website 
(NIST,  2003). 

Therefore  the  change  in  lithium  concentration  over 
time  can  be  related  to  the  change  in  neutron  intensity.  The 
negative  logarithm  of  the  ratio  of  two  measurements  of  the 
beam  intensity,  referred  to  as  the  optical  density  (OD),  is 
proportional  to  the  change  in  lithium  concentration 

OD  =  —In (I(t2)/I{h))  =  OuNA8{cLi(t2)-cLi{ti))  A cLi. 

(3) 

Note,  the  change  in  optical  density  can  only  be  related  to 
a  change  in  Lithium  concentration  in  the  solid  electrode, 
by  assuming  that  nothing  else  in  the  battery  is  changing  (or 
moving). 


Figure  3:  (Left)  Neutron  imaging  setup.  The  neutron  beam 
is  attenuated  as  it  passes  through  the  pouch  cell  to  reach  the 
detector,  producing  a  2-D  projection  of  the  battery  onto  the 
x-y  plane.  (Right)  Neutron  image  of  pouch  cell 


density  of  lithium  atoms  in  solution  is  too  low,  and  hence 
the  changes  are  imperceptible.  This  constraint  of  the  imag¬ 
ing  system  is  not  restrictive.  In-fact  this  constraint  prevents 
us  from  imaging  the  fast  changing  electrolyte  concentra¬ 
tion.  On  the  other  hand,  due  to  the  high  neutron  cross  sec¬ 
tion  of  the  hydrocarbon  based  solvents,  neutron  radiogra¬ 
phy  has  been  used  to  see  the  absorption  of  the  electrolyte 
into  the  battery  electrode  structure  during  cycling  (Lanz 
etal.,2001). 

Since  the  image  is  formed  by  neutron  counting,  which 
is  a  random  process,  there  is  some  variation  in  the  measure¬ 
ment.  A  reduction  of  the  noise  in  the  image  can  be  achieved 
by  increasing  the  exposure  time  of  the  image  and  counting 
more  neutrons,  but  then  the  decreased  temporal  resolution 
may  be  unsuitable  for  imaging  transient  behavior  such  as 
battery  charging.  Another  approach  to  reducing  the  mea¬ 
surement  variance  is  to  average  many  pixels  together  over 
a  region  with  uniform  characteristics,  in  this  case  lithium 
concentration  in  the  electrode.  By  treating  each  pixel  as  an 
independent  measurement  and  taking  a  spatial  average  we 
can  achieve  the  same  effect  as  increasing  the  exposure  time 
at  the  sacrifice  of  in-plane  resolution. 

The  planar  pouch-cell  construction  has  a  unique  ad¬ 
vantage  over  cylindrical  designs  for  neutron  imaging  as 
the  difference  between  anode  and  cathode  layers  are  eas¬ 
ily  identified  in  the  neutron  image  and  the  electrode  width 
8,  along  the  beam  path,  can  be  easily  tailored  to  optimize 
the  beam  transmission.  Since  we  are  interested  to  measure 
the  change  in  lithium  concentration  across  the  electrode 
(in  the  v-direction  or  through  plane),  and  the  concentration 
along  the  electrode  (along  the  y-axis  or  in-plane)  should 
be  uniform;  averaging  along  the  in-plane  direction  does 
not  sacrifice  any  information.  Therefore  the  use  of  line- 
profiles,  where  many  (approximately  1000)  pixels  along 
the  v-direction  are  averaged  together,  is  used  for  further 
data  analysis  and  discussion.  The  line  profile  shown  in 
Fig.  4,  is  calculated  based  on  1000  pixel  average  from  the 
actual  image  shown  in  Fig.  3  as  indicated  by  the  rectangular 
box. 


As  can  be  seen  in  Fig.  4,  there  is  considerable  variabil¬ 
ity  in  the  transmission  profile  of  each  folded  layer  in  the 
pouch  cell.  The  understanding  and  modeling  of  this  source 
of  variability  will  be  investigated  in  a  future  publication. 
In  the  following  paragraph  we  compare  the  first  two  folds 
with  a  modeled  transmission  line  profile  corresponding  to 
a  charged  battery.  Note  here  that  the  optical  density,  hence 
the  change  in  lithium  concentration  is  more  uniform  across 
the  multiple  folded  layers  despite  the  considerable  variabil¬ 
ity  in  the  transmission  line  profile  across  the  pouch  cell  as 
shown  later  in  Fig.  9. 

2.2  Image  Formation  and  Capture  Process 


Measurement  of  the  change  in  lithium  ion  concentra¬ 
tion  in  the  electrolyte  is  not  feasible  because  the  number 


The  thickness  of  commercial  electrode  materials  push 
the  resolution  limit  of  currently  available  high  resolution 


3 


UNCLASSIFIED 


pixel  number 


Figure  4:  Transmission  line-profile  for  the  reference 
(charged)  state,  calculated  by  averaging  over  the  region  of 
the  image  shown  in  Fig.  3.  Each  pixel  is  5  /im. 

cross-strip  anode  Micro-Channel  Plate  (MCP)  detectors 
(Tremsin  et  al.,  2005;  Siegmund  et  al.,  2009;  Hussey  et  al., 
2010).  Prior  attempts  to  investigate  Lithium  ion  batteries 
using  neutron  imaging  (Kamata  et  al.,  1997)  were  limited 
by  the  resolution  of  the  imaging  system,  but  the  new  high 
resolution  MCP  detector,  with  10  lim  pore  spacing,  is  capa¬ 
ble  of  resolving  features  on  the  order  of  1 3.5/im.  The  MCP 
has  6/im  diameter  pores  on  1 0/i m  centers,  with  a  5, lim 
pixel  pitch.  This  high  resolution  detector  enables  the  mea¬ 
surement  of  lithium  concentration  across  the  solid  phase  of 
the  battery  electrodes,  as  the  typical  electrode  thickness  is 
between  40-100/im  (Doyle  et  al.,  1993;  Smith  et  al.,  2010; 
Speltino  et  al.,  2009). 

Figure  5  shows  a  schematic  of  a  lithium-ion  pouch  cell 
and  corresponding  simulation  of  the  resulting  ideal  trans¬ 
mission  line  profile  according  to  the  Lambert-Beer  law 
Eq.  (1).  The  cell  consists  of  alternating  double  sided  cur¬ 
rent  collectors;  LiFePCU  on  aluminum  for  the  positive  and 
carbon  on  copper  for  the  negative  electrodes  as  shown  in 
the  unwrapped  pouch  cell  construction  of  Fig.  1  and  in 
the  transmission  line  profiles  with  the  lower  schematic  of 
Fig.  5.  The  electrodes  are  separated  by  a  porous  polyethy¬ 
lene  or  polypropylene  separator  material.  The  multi-layer 
sheets  of  electrode  material  are  then  folded  10  times  to  pro¬ 
duce  a  10  layer  pouch  cell  shown  in  the  neutron  image  of 
Fig.  3  and  line  profile  in  Fig.  4.  The  peaks  in  both  the  ideal 
and  blurred  line  profiles  are  caused  by  the  aluminum  foil 
current  collector  which  has  an  order  of  magnitude  lower 
neutron  attenuation  than  other  materials  in  the  pouch  cell. 
The  highest  attenuation  is  expected  to  come  from  the  sep¬ 
arator  region,  since  the  membrane  consists  of  plastic  mate¬ 
rial  that  strongly  attenuates  neutrons,  however  it  is  not  easy 
to  identify  in  either  the  images  or  line-profile  due  to  the 
thinness  of  the  membrane  (  ss  25/im)  and  proximity  to  the 


active  material. 

The  line  profile  associated  to  the  fully  charged  state 
and  the  material  and  electrode  thickness  is  shown  in  Fig.  5 
with  the  dash-dot  line.  Using  modeling  results  from 
(Domenico  et  al.,  2010)  a  uniform  Li  concentration  across 
the  electrode  is  expected  due  to  the  low  C-rate.  The 
negative  electrode  thickness  is  approx  37  /im  (each  side) 
graphite  with  10  /im  copper  foil  backing,  the  positive  elec¬ 
trode  is  53  ji m  (each  side)  with  20  fim  aluminum  foil  back¬ 
ing,  the  separators  are  25  /im  (PP  or  PE).  The  peaks  in  the 
both  the  ideal  and  blurred  line  profiles  correspond  the  alu¬ 
minum  foil  current  collector  which  has  an  order  of  magni¬ 
tude  lower  neutron  attenuation  than  other  materials  in  the 
pouch  cell.  Proprietary  information  on  the  material  con¬ 
struction  and  composition  was  used  to  model  the  ideal  neu¬ 
tron  transmission  through  each  layer. 

The  image  formation  process  introduces  several  arti¬ 
facts  changing  the  ideal  line  profile  predicted  by  Eq.  (1) 
shown  by  the  dash-dot  line  into  the  smoother  line  profile 
shown  by  the  solid  and  the  dashed  lines.  This  smooth¬ 
ing  process  complicates  the  quantification  of  the  lithium 
concentration  from  the  detected  signal.  These  artifacts  are 
caused  by  4  things;  improper  alignment  of  the  battery  with 
respect  to  the  neutron  beam,  geometric  unsharpness  due  to 
the  pinhole  imaging  system,  roughness  or  curvature  of  the 
object  (non-parallel  layers),  and  blurring  of  the  signal  in¬ 
troduced  by  the  detector.  Proper  alignment  of  the  battery 
layers,  parallel  to  the  neutron  beam  direction,  is  critical  to 
produce  the  expected  transmission  line  profile.  If  the  bat¬ 
tery  is  slightly  misaligned  (rotation  around  the  y-axis)  the 
transmission  line  profile  corresponding  to  the  various  lay¬ 
ers  in  the  battery  may  appear  blurred,  as  shown  by  the  line 
profile  corresponding  to  the  rotated  image  in  Fig.  5,  due 
to  neutron  transmission  through  two  or  more  layers  of  the 
pouch  cell. 

The  expected  line  profile  is  then  convolved  with  a  point 
spread  function  (PSF),  to  simulate  the  blurring  of  the  detec¬ 
tor,  and  imaging  optics  (Hussey  et  al.,  2010).  Deconvolu¬ 
tion  of  the  data,  to  reconstruct  the  ideal  line  profile  is  dif¬ 
ficult  as  this  would  intensify  the  noise  present  in  the  signal 
(Dupe  et  al.,  2008),  yielding  a  meaningless  result.  There¬ 
fore  we  simulate  the  blurring  caused  by  the  measurement 
system  and  filter  the  modeled  transmission  line  profile,  in 
order  to  compare  the  measured  data  with  the  expected  re¬ 
sults. 

The  geometry  of  the  imaging  setup,  shown  in  Fig.  3, 
contributes  to  blurring  of  the  resulting  signal.  This  geo¬ 
metric  unsharpness  can  be  modeled  by  convolving  the  ideal 
line  profile  with  a  rectangular  window  of  width  I/?  =  zD/L , 
where  L/D  =  3000  is  the  ratio  of  source  to  detector  dis¬ 
tance  over  the  aperture  diameter,  and  z  is  the  distance  be¬ 
tween  the  detector  and  the  object  which  is  imaged  (Hussey 
et  al.,  2010;  A. A.  Harms,  1986).  For  our  experimental 
setup  z  =  5.4  cm,  so  the  theoretical  geometric  unsharpness 
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Figure  5:  A  model  of  the  transmission  line  profiles  base 
on  the  layer  structure  of  the  Li-ion  pouch  cell  (bottom)  is 
shown  with  dash-dot  corresponding  to  the  ideal  beam  in¬ 
tensity,  solid  to  the  modeled  transmission  line  profile  after 
a  0.09  °  rotation  of  the  pouch  cell  with  respect  to  the  beam, 
and  dashed  to  the  model  of  blurring  associated  with  neutron 
capture  process  in  the  detector. 


is  L/g  =  18  11  m.  The  PSF  corresponding  to  the  detector  can 
be  modeled  by  a  Gaussian  kernel  of  width  Od  =  9.2  /im, 
taken  from  measurements  of  the  detector  (Hussey  et  ah, 
2010).  The  resulting  simulated  line  profile,  which  accounts 
for  the  three  sources  of  image  degradation,  is  shown  in 
Fig.  5.  The  simulated  line  profile  with  black  dashed  line 
can  be  compared  with  the  actual  (measured)  line  profile  in 
Fig.  4. 

3.  NEUTRON  IMAGING  EXPERIMENTAL 
RESULTS 

This  section  describes  the  testing  protocol  used  to 
stimulate  the  battery  and  the  steady  state  conditions  under 
which  neutron  images  were  acquired.  Low  current  rates 
were  applied  to  the  battery  in  order  to  ensure  a  flat  distri¬ 
bution  of  lithium  concentration  across  the  solid  electrode 
(Domenico  et  ah,  2010),  which  could  be  easily  measured 
and  facilitate  comparison  with  the  theoretical  results. 

3.1  Experimental  Procedure 

The  lower  subplot  of  Fig.  6  shows  the  current  profile 
applied  to  the  battery  (dashed  line),  and  the  stored  charge 
in  mAh  (solid  line).  The  resulting  voltage  profile  is  shown 
in  the  upper  subplot.  The  initial  state  of  the  battery  was 
discharged,  and  at  rest.  The  battery  was  aligned  and  im¬ 
aged  before  charging  at  C/5  rate  for  2.5  hours  to  50%  SOC. 
After  another  set  of  images  was  acquired,  the  battery  was 


charged  to  100  %  SOC  using  a  constant-current,  constant- 
voltage  (3.6V)  charging  algorithm.  Neutron  images  of  the 
battery  in  a  known  reference  (fully  charged)  state  were  then 
acquired;  where  most  of  the  mobile  lithium  ions  are  ex¬ 
pected  to  be  intercalated  into  the  carbon  or  negative  elec¬ 
trode.  Imaging  concluded  at  t  =  940  (min)  indicated  in 
Fig.  6.  Several  image  frames  are  combined  (in  time)  to 
produce  an  image  with  a  very  long  effective  exposure  time 
(1  hour).  Longer  exposure  times  are  utilized  to  reduce  the 
uncertainty  in  the  reference  image  from  neutron  counting 
statistics. 


Figure  6:  Charge  and  discharge  profile  applied  to  the  bat¬ 
tery  during  imaging. 

The  battery  was  then  discharged  to  50%  SOC  at  C/5 
rate  and  imaged  during  a  period  of  rest  for  1  hour  with  zero 
applied  current,  at  time  t$  indicated  in  Fig.  6.  Once  the 
imaging  was  completed,  the  discharge  continued  to  3V,  or 
0%  SOC,  and  a  final  set  of  images  was  acquired  under  rest 
conditions  at  time  t&. 

3.2  Expansion  of  electrode  materials  (swelling  during 
charge) 

Thermal  expansion  of  the  materials,  and/or  expansion 
of  anode  solid  electrode  material  during  lithium  intercala¬ 
tion  (Fuller  et  ah,  1994),  requires  special  attention  to  cor¬ 
rectly  interpret  the  neutron  radiographs.  The  change  in 
lithium  concentration  is  proportional  to  the  negative  log¬ 
arithm  of  the  ratio  of  two  images  (Hussey  et  ah,  2010), 
assuming  nothing  else  within  the  image  is  moving.  Move¬ 
ment  of  the  active  material,  as  suggested  by  the  peak  mis¬ 
alignment  shown  in  Fig.  7  may  lead  to  erroneous  estimate 
of  the  change  in  lithium  concentration.  Therefore  methods 
developed  in  subsection  3.3  are  used  to  estimate  and  cor¬ 
rect  for  the  expansion/contraction  of  the  pouch  cell  during 
charging/discharging. 

The  pouch  cell  is  secured  to  an  aluminum  fixture  on  the 
right  side  of  the  image,  and  therefore  material  expansion 
will  produce  a  shift  with  increasing  displacement  (from 
right  to  left  or  decreasing  pixel  number)  as  the  layers  each 
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Figure  7:  The  two  left  most  peaks  from  the  line  profile  of 
neutron  transmission  through  the  battery,  from  Fig.  4,  are 
shown  with  a  schematic  of  the  pouch  cell  showing  only 
two  layers.  Expansion  of  the  active  material  during  charg¬ 
ing  is  visible  in  the  peak  shift,  c.f.  the  dash-dot  and  dashed 
lines,  corresponding  to  the  left-most  peak,  which  is  furthest 
from  the  mounting  plate  and  therefore  subject  to  the  great¬ 
est  translational  motion.  The  total  thickness  of  each  unit 
battery  is  126+84+50=260  jtim,  which  agrees  well  with  the 
average  peak  spacing  of  52.5  pixels  or  262.5  jtim 


move  to  the  left,  and  expand.  The  measured  peak  displace¬ 
ments  are  shown  in  Fig.  8  and  indeed  follow  the  linearly 
increasing  shift  pattern  with  the  distance  away  from  the 
fixed  mount.  The  average  displacement  corresponds  to  a 
3%  length  change  if  attributed  solely  to  the  anode  active 
material,  which  is  within  the  range  of  0-7%  predicted  by 
(Fuller  et  ah,  1994). 


Figure  8:  Peak  displacements  for  each  of  the  10  layers, 
from  the  fully  charged  reference  state  to  the  discharged 
state.  Shrinkage  of  the  active  material  is  observed  upon 
discharging.  Each  pixel  corresponds  to  5  jum. 


3.3  Corrected  Line  Profiles  and  Optical  Density 

The  peaks  in  the  line  profile,  which  correspond  to 
the  aluminum  current  collector  for  the  positive  electrode, 
provide  the  most  well  defined  feature  for  determining  the 
amount  of  translational  motion  present  in  the  image.  Due 
to  the  noise  present  in  the  image  and  in  the  line  profile,  es¬ 
timation  of  the  peak  height  and  location  is  necessary,  rather 
than  simply  finding  the  maximum  value.  An  estimate  of  the 
peak  location  can  be  found  by  fitting  a  Gaussian  kernel. 
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locally  to  the  data  using  least-squares  algorithm  with  four 
parameters  (cr,  x,  a ,  b )  to  determine  the  peak  location,  x, 
with  sub-pixel  resolution.  The  reference  line  profile  is  then 
stretched,  using  an  affine  transformation  (linear  interpola¬ 
tion)  to  match  the  peak  locations.  This  linear  scaling  of  the 
reference  transmission  line  profile  is  necessary  for  calculat¬ 
ing  the  optical  density  using  Eq.  (3),  shown  in  Fig.  9,  so  that 
the  same  regions  of  the  active  material  are  compared  and 
valid  results  can  be  obtained.  Despite  the  non-uniformity 
in  the  detected  transmission  line  profile  shown  in  Fig.  4, 
the  optical  density  from  Eq.  (3)  shows  fairly  uniform  char¬ 
acteristics  between  pouch  layers.  A  zoomed  plot  showing 
the  optical  density  first  two  layers  is  shown  in  Fig.  10.  In 
this  figure  the  optical  density  of  different  material  layers 
within  the  unit  cell  are  not  clearly  visible.  The  optical  den¬ 
sity  is  proportional  to  lithium  concentration;  thus  the  only 
expected  change  in  optical  density  should  occur  in  the  re¬ 
gions  corresponding  to  the  anode  and  cathode  active  mate¬ 
rial.  The  reference  state  is  taken  to  be  fully  charged,  so  that 
lithium  accumulation  (positive  values  of  optical  density) 
occurs  in  the  cathode  layer  when  the  battery  is  discharged 
relative  to  the  reference,  and  the  negative  electrode  is  los¬ 
ing  lithium  therefore  the  optical  density  decreases.  The  two 
solid  lines  in  the  optical  density  plot  corresponding  to  50% 
SOC  are  not  identical.  This  difference,  however  can  be  ex¬ 
plained  because  the  first  set  of  images  (corresponding  to 
the  battery  charging  cycle)  were  taken  before  a  change  in 
the  nuclear  reactor  which  slightly  changes  the  optical  axis 
of  the  neutron  beam. 


3.4  Theoretical  vs  Experimental  Results 

This  section  shows  a  comparison  between  the  modeled 
optical  density  corresponding  to  a  change  in  lithium  con¬ 
centration  predicted  by  the  battery  model  (Speltino  et  ah, 
2009;  Domenico  et  ah,  2010)  with  the  experimental  results. 
Despite  their  complexity,  battery  models  based  on  electro¬ 
chemistry  laws  such  as  (Speltino  et  ah,  2009;  Domenico 
et  ah,  2010;  Fuller  et  ah,  1994;  Santhanagopalan  et  ah, 
2008)  are  generally  preferred  to  the  equivalent  circuit  or 
to  other  kinds  of  simplified  models,  because  they  also  pre¬ 
dict  physical  battery  states,  which  have  a  relevant  effect  in 
the  automotive  application,  where  the  battery  often  expe¬ 
riences  the  stress  of  very  high  transient  loads  as  discussed 
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pixel  number 

Figure  9:  Line  profile  of  the  optical  density,  which  is 
proportional  to  lithium  concentration,  at  various  states  of 
charge.  The  reference  state  is  taken  to  be  fully  charged, 
so  that  lithium  accumulation  (positive  values)  occurs  in  the 
cathode  layer  when  the  battery  is  discharged  relative  to  the 
reference. 
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Figure  10:  Zoomed  plot  showing  the  optical  density  cal¬ 
culated  from  radiographic  measurement  using  a  1000  px 
averaged  line  profile,  corresponding  to  the  first  two  battery 
layers. 

in  (Smith  and  Wang,  2006).  In  order  to  simplify  the  prob¬ 
lem  and  determine  the  quality  of  the  metrology  technique 
the  battery  is  imaged  during  rest  conditions.  The  electro¬ 
chemical  model  is  then  used  to  predict  to  the  steady  state 
lithium  distribution  in  the  solid  electrode  cse(x),  which  is  a 
constant  flat  profile  across  the  electrode. 


Figure  11:  Ideal  simulated  transmission  line  profile, 
blurred  line  profile  and  resulting  optical  densities  in  the  dis¬ 
charged  state,  with  the  reference  taken  at  the  fully  charged 
state. 

The  experimental  results  show  reasonable  matching 
with  the  numerical  simulation  for  the  optical  density  ob¬ 
tained  from  the  simulation  of  the  lithium  distribution  in 
the  solid  electrodes  during  cycling  of  the  battery  at  low  C- 
rate  and  imaging  during  steady  state.  The  ideal  line  pro¬ 
file,  corresponding  to  the  change  in  lithium  concentration 
calculated  directly  from  the  line  profile  using  the  Lambert  - 
Beer  law,  does  not  agree  with  the  measured  data;  however, 
accounting  for  the  artifacts  of  the  imaging  system  show 
much  better  agreement  with  experimental  results.  A  higher 
value  of  <Jd  =  30  microns,  than  the  theoretical  value  of 
0/3  =9.2  microns,  was  needed  to  match  the  experimen¬ 
tal  data.  The  higher  0/>  causes  increased  blurring  of  the 
simulated  line  profile  and  shows  better  agreement  with  the 
blurred  peaks  in  the  transmission  line  profile,  shown  in 
Fig.  1 1;  in  particular  the  peak  near  pixel  #38,  correspond¬ 
ing  to  the  anode  active  material  in  the  discharged  state, 
which  becomes  flat  during  charged  conditions.  The  in¬ 
creased  blurring  observed  in  the  neutron  images  may  be 
attributed  to  in-plane  roughness  or  curvature  of  the  bat¬ 
tery  layers,  which  highlights  the  importance  of  designing 
battery  cells  which  are  suitable  for  imaging.  Quantitative 
assessments  of  the  change  in  lithium  concentration  distri¬ 
bution  across  the  electrode  requires  calibration  of  the  mea¬ 
surement  tool  using  a  Lithium  Fluoride  standard  by  NIST 
neutron  imaging  facility,  which  is  currently  in  progress  and 
once  completed  will  allow  the  translation  of  measured  op¬ 
tical  density  into  the  change  of  lithium  concentration  in  the 
solid  electrode. 

CONCLUSIONS 

The  experimental  results  show  reasonable  matching 
with  numerical  simulation  for  the  optical  density  obtained 
from  a  model  of  the  lithium  distribution  in  the  solid  elec¬ 
trodes  during  cycling  of  the  battery  at  low  C-rate.  Attempts 
to  calculate  the  change  in  lithium  concentration  directly 
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from  the  line  profile  using  the  Lambert-Beer  law  demon¬ 
strated  the  need  to  consider  artifacts  of  the  imaging  system 
when  trying  to  make  quantitative  assessments  of  the  change 
in  lithium  concentration  distribution  across  the  electrode. 
The  concentration  measurements  can  then  be  used  to  val¬ 
idate  models  for  battery  management  in  applications  that 
require  high  power  demand  and  hence  need  accurate  esti¬ 
mates  of  the  both  the  bulk  and  the  spatiotemporal  Li  con¬ 
centration. 
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